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zr> Multiphase Coriolis Flowmeter 

^ CROSS-REFERENCE TO RELATED APPLICATIONS 

This application claims priority to U.S. Application Serial No. 60/445,795, filed on 
February 10, 2003, and titled MULTIPHASE CORIOLIS FLOWMETER. This application also 
claims priority to U.S. Application Serial No. 60/452,934, filed on March 10, 2003, and titled 
5 MULTIPHASE CORIOLIS FLOWMETER. Both of these applications are hereby incorporated 
by reference. 



TECHNICAL FIELD 

This description relates to flowmeters. 



BACKGROUND 

10 Flowmeters provide information about materials being transferred through a conduit, or 

flowtube. For example, mass flowmeters provide an indication of the mass of material being 
transferred through a conduit. Similarly, density flowmeters, or densitometers, provide an 
indication of the density of material flowing through a conduit. Mass flowmeters also may 
provide an indication of the density of the material, and therefore an indication of the volumetric 

15 flow rate. 

For example, Coriolis-type mass flowmeters are based on the Coriolis effect, in which 
material flowing through a conduit becomes a radially-travelling mass that is affected by a 
Coriolis force and therefore experiences an acceleration. Many Coriolis-type mass flowmeters 
induce a Coriolis force by sinusoidally oscillating a conduit about a pivot axis orthogonal to the 
20 length of the conduit. In such mass flowmeters, the Coriolis reaction force experienced by the 
traveling fluid mass is transferred to the conduit itself and is manifested as a deflection or offset 
of the conduit in the direction of the Coriolis force vector in the plane of rotation. 

SUMMARY 

According to one general aspect, a flowmeter includes a vibratable flowtube, a driver 
25 connected to the flowtube and operable to impart motion to the flowtube, a sensor connected to 
the flowtube and operable to sense the motion of the flowtube and generate a sensor signal, and a 
controller connected to receive the sensor signal, the controller being operable to determine a 

1 



Attorney Docket No. 1 2780-02700 1 



first flow rate of a first phase within a two-phase flow through the flowtube and determine a 
second flow rate of a second phase within the two-phase flow. 

Implementations may include one or more of the following features. For example, the 
first phase may include a gas and the second phase may include a liquid. 
5 The controller may be operable to input an apparent density of the two-phase flow 

detected by the flowmeter and output a corrected density of the two-phase flow. The controller 
may be operable to correct the apparent density based on a theoretical relationship between the 
apparent density and the corrected density, or based on an empirical relationship between the 
apparent density and the corrected density (such as, for example, a table storing relationships 
10 between the apparent density and the corrected density). 

The controller may be operable to input an apparent mass flow rate of the two-phase flow 
detected by the flowmeter and output a corrected mass flow rate of the two-phase flow. The 
controller may be operable to correct the apparent mass flow rate based on a theoretical or 
empirical relationship, such as a tabular relationship, between the apparent mass flow rate and 
15 the corrected mass flow rate. 

The controller may be operable to input an apparent first phase fraction of the two-phase 
flow detected by the flowmeter that defines an amount of the first phase in the two-phase flow 
and output a corrected first phase fraction of the two-phase flow. The controller may be operable 
to input a first phase fraction of the two-phase flow detected by a phase fraction sensor that is 
20 external to the flowmeter. 

The controller may be operable to determine the first flow rate and the second flow rate 
based on corrected values for a detected density and detected mass flow rate of the two-phase 
flow. The controller may be operable to determine the first flow rate and the second flow rate 
based on a corrected value for a detected first phase fraction that defines an amount of the first 
25 phase in the two-phase flow. The controller may be operable to determine the first flow rate and 
the second flow rate based on densities of the first phase and the second phase, respectively. 

The controller may be operable to determine a first superficial velocity of the first phase 
and a second superficial velocity of the second phase, based on the first flow rate and the second 
flow rate, respectively. The controller may be operable to determine a flow regime of the two- 
30 phase flow, based on the first superficial velocity and the second superficial velocity. The 

controller may be operable to determine a slip velocity between the first phase and the second 
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phase, based on an average velocity of the first phase and an average velocity of the second 
phase. The controller may be operable to provide corrections to the first flow rate and the second 
flow rate, based on the first and second superficial velocities, the determined flow regime, or the 
slip velocity, to thereby obtain a corrected first flow rate and a corrected second flow rate. 
5 According to another general aspect, a method includes determining a bulk density of a 

two-phase flow through a flowtube, the two-phase flow including a first phase and a second 
phase, determining a bulk mass flow rate of the two-phase flow, and determining a first mass 
flow rate of the first phase, based on the bulk density and the bulk mass flow rate. 

Implementations may include one or more of the following features. For example, a 
10 second mass flow rate of the second phase may be determined, based on the bulk density and the 
bulk mass flow rate. In determining the bulk density, an apparent bulk density of the two-phase 
flow may be determined, and the apparent bulk density may be corrected to obtain the bulk 
density. 

In correcting the apparent bulk density, the apparent bulk density may be input into a 
15 theoretical relationship that relates the apparent bulk density to a corrected bulk density, or may 
be input into an empirical relationship that relates the apparent bulk density to a corrected bulk 
density. 

In correcting the apparent bulk density, a first density of the first phase may be input. A 
first phase fraction of the two-phase flow may be determined, based on the bulk density, the first 
20 density of the first phase, and a second density of the second phase. In determining the first mass 
flow rate of the first phase, the first mass flow rate may be determined based on the first phase 
fraction and the first density. 

A first superficial velocity of the first phase and a second superficial velocity of the 
second phase may be determined, based on the first mass flow rate and the second mass flow 
25 rate, respectively. A flow regime of the two-phase flow may be determined, based on the first 
superficial velocity and the second superficial velocity. A slip velocity between the first phase 
and the second phase may be determined, based on an average velocity of the first phase and an 
average velocity of the second phase. Corrections may be provided to the first flow rate and the 
second flow rate, based on the first and second superficial velocities, the determined flow 
30 regime, or the slip velocity. 

The first phase may include a gas and the second phase may include a liquid. 
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According to another general aspect, a flowmeter controller includes a density correction 
system operable to input an apparent density of a two-phase flow and output a corrected density 
of the two-phase flow, the two-phase flow including a first phase and a second phase, a mass 
flow rate correction system operable to input an apparent mass flow rate of the two-phase flow 
5 and output a corrected mass flow rate of the two-phase flow, and a flow component mass flow 
rate determination system operable to determine a first mass flow rate of the first phase, based on 
the corrected density and the corrected mass flow rate. 

Implementations may include one or more of the following features. For example, the 
flow component mass flow rate determination system may be operable to determine a second 
10 mass flow rate of the second phase, based on the corrected density and the corrected mass flow. 

The first phase may include a liquid and the second phase may include a gas. A phase 
fraction determination system may be included that is operable to determine a corrected phase 
fraction of the two-phase flow, wherein the flow component mass flow rate determination system 
may be operable to determine the first flow rate and the second flow rate based on the corrected 
15 phase fraction. The phase fraction determination system may be a void fraction determination 
system that determines an amount of the gas in the two-phase flow. 

A superficial velocity determination system may be included that is operable to determine 
a first superficial velocity of the first phase and a second superficial velocity of the second phase. 
The flowmeter controller may include a flow regime determination system operable to determine 
20 a flow regime of the two-phase flow. 

The flow regime determination system may be further operable to determine a phase slip 
velocity with respect to an average velocity of the first phase and an average velocity of the 
second phase. The flow component mass flow rate determination system may be operable to 
improve the determination of the first mass flow rate and the second mass flow rate, based on the 
25 first and second superficial velocities, the flow regime, or the phase slip velocity. 

The details of one or more implementations are set forth in the accompanying drawings 
and the description below. Other features will be apparent from the description and drawings, 
and from the claims. 
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DESCRIPTION OF DRAWINGS 

FIG 1 A is an illustration of a Coriolis flowmeter using a bent flowtube. 
FIG IB is an illustration of a Coriolis flowmeter using a straight flowtube. 
FIG 2 is a block diagram of a Coriolis flowmeter. 

FIG 3 is a flowchart illustrating an operation of the coriolis flowmeter of FIG 2. 
FIG 4 is a flowchart illustrating techniques for determining liquid and gas flow rates for a 
two-phase flow. 

FIGS. 5A and 5B are graphs illustrating a percent error in a measurement of void fraction 
and liquid fraction, respectively. 

FIG 6 is a graph illustrating a mass flow error as a function of a drop in density for a 
flowtube having a particular orientation and over a selected flow range. 

FIG 7 is a flowchart illustrating techniques for correcting density measurements. 

FIG 8 is a table showing a relationship between an apparent density drop and an apparent 
mass flow rate of the two-phase flow. 

FIG 9 is a flowchart illustrating techniques for determining void fraction measurements. 

FIG 10 is a flowchart illustrating techniques for determining corrected mass flow rate 
measurements. 

FIG 11 is a table showing a relationship between an apparent mass flow rate and a 
corrected density drop of the two-phase flow. 

FIGS. 12-14 are graphs illustrating examples of density corrections for a number of 
flowtubes. 

FIGS. 15-20 are graphs illustrating examples of mass flow rate corrections for a number 
of flowtubes. 

DETAILED DESCRIPTION 

Types of flowmeters include digital flowmeters. For example, U.S. Patent 6,31 1,136, 
which is hereby incorporated by reference, discloses the use of a digital flowmeter and related 
technology including signal processing and measurement techniques. Such digital flowmeters 
may be very precise in their measurements, with little or negligible noise, and may be capable of 
enabling a wide range of positive and negative gains at the driver circuitry for driving the 
conduit. Such digital flowmeters are thus advantageous in a variety of settings. For example, 
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commonly-assigned U.S. Patent 6,505,519, which is incorporated by reference, discloses the use 
of a wide gain range, and/or the use of negative gain, to prevent stalling and to more accurately 
exercise control of the flowtube, even during difficult conditions such as two-phase flow (e.g., a 
flow containing a mixture of liquid and gas). 
5 Although digital flowmeters are specifically discussed below with respect to, for 

example, FIGS. 1 and 2, it should be understood that analog flowmeters also exist. Although 
such analog flowmeters may be prone to typical shortcomings of analog circuitry, e.g., low 
precision and high noise measurements relative to digital flowmeters, they also may be 
compatible with the various techniques and implementations discussed herein. Thus, in the 

10 following discussion, the term "flowmeter" or "meter" is used to refer to any type of device 
and/or system in which a Coriolis flowmeter system uses various control systems and related 
elements to measure a mass flow, density, and/or other parameters of a material(s) moving 
through a flowtube or other conduit. 

FIG. 1 A is an illustration of a digital flowmeter using a bent flowtube 102. Specifically, 

15 the bent flowtube 102 may be used to measure one or more physical characteristics of, for 
example, a (traveling) fluid, as referred to above. In FIG. 1 A, a digital transmitter 104 
exchanges sensor and drive signals with the bent flowtube 102, so as to both sense an oscillation 
of the bent flowtube 102, and to drive the oscillation of the bent flowtube 102 accordingly. By 
quickly and accurately determining the sensor and drive signals, the digital transmitter 104, as 

20 referred to above, provides for fast and accurate operation of the bent flowtube 102. Examples 
of the digital transmitter 104 being used with a bent flowtube are provided in, for example, 
commonly-assigned U.S. Patent 6,31 1,136. 

FIG. IB is an illustration of a digital flowmeter using a straight flowtube 106. More 
specifically, in FIG. IB, the straight flowtube 106 interacts with the digital transmitter 104. Such 

25 a straight flowtube operates similarly to the bent flowtube 102 on a conceptual level, and has 
various advantages/disadvantages relative to the bent flowtube 102. For example, the straight 
flowtube 106 may be easier to (completely) fill and empty than the bent flowtube 102, simply 
due to the geometry of its construction. In operation, the bent flowtube 102 may operate at a 
frequency of, for example, 50-1 10 Hz, while the straight flowtube 106 may operate at a 

30 frequency of, for example, 300-l,000Hz. The bent flowtube 102 represents flowtubes having a 
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variety of diameters, and may be operated in multiple orientations, such as, for example, in a 

vertical or horizontal orientation. 

Referring to FIG. 2, a digital mass flowmeter 200 includes the digital transmitter 104, one 

or more motion sensors 205, one or more drivers 210, a flowtube 215 (which also may be 
5 referred to as a conduit, and which may represent either the bent flowtube 102, the straight 

flowtube 106, or some other type of flowtube), and a temperature sensor 220. The digital 

transmitter 104 may be implemented using one or more of, for example, a processor, a Digital 

Signal Processor (DSP), a field-programmable gate array (FPGA), an ASIC, other programmable 

logic or gate arrays, or programmable logic with a processor core. It should be understood that, 
10 as described in 6,31 1,136, associated digital-to-analog converters may be included for operation 

of the drivers 210, while analog-to-digital converters may be used to convert sensor signals from 

the sensors 205 for use by the digital transmitter 104. 

The digital transmitter 104 generates a measurement of, for example, density and/or mass 

flow of a material flowing through the flowtube 215, based at least on signals received from the 
15 motion sensors 205. The digital transmitter 104 also controls the drivers 210 to induce motion in 

the flowtube 215. This motion is sensed by the motion sensors 205. 

Density measurements of the material flowing through the flowtube are related to, for 

example, the frequency of the motion of the flowtube 215 that is induced in the flowtube 215 by 

a driving force supplied by the drivers 210, and/or to the temperature of the flowtube 215. 
20 Similarly, mass flow through the flowtube 215 is related to the phase and frequency of the 

motion of the flowtube 215, as well as to the temperature of the flowtube 215. 

The temperature in the flowtube 215, which is measured using the temperature sensor 

220, affects certain properties of the flowtube, such as its stiffness and dimensions. The digital 

transmitter 104 may compensate for these temperature effects. Also in FIG. 2, a pressure sensor 
25 225 is in communication with the transmitter 104, and is connected to the flowtube 215 so as to 

be operable to sense a pressure of a material flowing through the flowtube 215. 

It should be understood that both the pressure of the fluid entering the flowtube 215 and 

the pressure drop across relevant points on the flowtube may be indicators of certain flow 

conditions. Also, while external temperature sensors may be used to measure the fluid 
30 temperature, such sensors may be used in addition to an internal flowmeter sensor designed to 

measure a representative temperature for flowtube calibrations. Also, some flowtubes use 
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multiple temperature sensors for the purpose of correcting measurements for an effect of 
differential temperature between the process fluid and the environment (e.g., a case temperature 
of a housing of the flowtube). As discussed in more detail below, one potential use for the inlet 
fluid temperature and pressure measurements is to calculate the actual densities of a liquid and 
5 gas in a two-phase flow, based on predefined formulae. 

A liquid fraction probe 230 refers to a device for measuring a volume fraction of liquid, 
e.g., water, when a liquid in the flowtube 215 includes water and another fluid, such as oil. Of 
course, such a probe, or similar probes, may be used to measure the volume fraction of a fluid 
other than water, if such a measurement is preferred or if the liquid does not include water. In 

10 the below description, a measured liquid is generally assumed to be water for the purposes of 
example, so that the liquid fraction probe 230 is generally referred to as a water fraction probe 
230, or a water-cut probe 230. 

A void fraction sensor 235 measures a percentage of a material in the flowtube 215 that is 
in gaseous form. For example, water flowing through the flowtube 215 may contain air, perhaps 

15 in the form of bubbles. Such a condition, in which the material flowing through the flowtube 
215 contains more than one material is generally referred to as "two-phase flow." In particular, 
the term "two-phase flow" may refer to a liquid and a gas; however, "two-phase flow" also may 
refer to other combinations of materials, such as two liquids (e.g., oil and water). 

Various techniques, represented generally in FIG. 2 by the void fraction sensor 235, exist 

20 for measuring the gas void fraction in a two-phase flow of liquid and gas. For example, various 
sensors or probes exist that may be inserted into the flow to determine a gas void fraction. As 
another example, a venturi tube (i.e., a tube with a constricted throat that determines fluid 
pressures and velocities by measurement of differential pressures generated at the throat as a 
fluid traverses the tube), relying on the fact that gas generally moves with a higher velocity than 

25 liquid(s) through a restriction, may be used to determine a pressure gradient and thereby allow a 
determination of the gas void fraction. Measurements of gas void fractions also may be obtained 
using equipment that is wholly external to the flowtube. For example, sonar measurements may 
be taken to determine gas void fraction. As a specific example of such a sonar-based system, the 
SONARtracTM gas void fraction monitoring system produced by CiDRA Corporation of 

30 Wallingford, Connecticut may be used. 
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In this description, an amount of gas in a flowing fluid, measured by the void fraction 
sensor or otherwise determined, is referred to as void fraction or a, and is defined as a = volume 
of gas / total volume = volume of gas / ( volume of liquid + volume of gas). Accordingly, a 
quantity referred to herein as the liquid fraction is defined as 1- a. 
5 In many applications where mass flow measurements are required, the void fraction of 

the flow can be as high as 20, 30, 40% or more. However, even at very small void fractions of 
0.5%, the fundamental theory behind the coriolis flowmeter becomes less applicable. 

Moreover, a presence of gas in the fluid flow also may affect a measurement of a density 
of the fluid flow, generally causing the density measurement to read lower. That is, it should be 

10 understood that a density pij qu id of a liquid flowing by itself through a flowtube will be higher 

than an actual density p true of a two-phase flow containing the liquid and a gas, since a density of 
the gas (e.g., air) will generally be lower than a density of the liquid (e.g., water) in the two- 
phase flow. In other words, there is a density reduction when gas is added to a liquid flow that 
previously contained only the liquid. 

15 Beyond this physical phenomenon, a coriolis meter measuring a two-phase fluid flow 

containing gas may output a density reading p app arent that is an ostensible measurement of the 
bulk density of the two-phase flow (e.g., of the water and air combined). This raw measurement 
Papparent will generally be different (lower) than the actual bulk density p^e of the two-phase 
flow. For example, the resonant frequency used by the flowmeter may be artificially high, due to 

20 relative motion of the gas in the fluid flow, which would cause the density measurement to read 
low. It should be understood that many conventional prior art flowmeters were unconcerned 
with this problem, since most such coriolis meters fail to continue operating (e.g. stall or output 
inaccurate measurements) at even the slightest amounts of void fraction. 

U.S. Patent No. 6,505,519, which is incorporated by reference above, discloses that such 

25 a variation of p app arent (i.e., an indicated raw or bulk density reading of a two-phase flow that is 
output by a coriolis flowmeter) from p^e (i.e., an actual raw or bulk density of the two-phase 
flow) may be characterized by a variety of techniques. As a result, a measured p ap parent may be 
corrected to obtain an actual bulk density p CO rrected, which is, at least approximately, equal to p^e- 
Somewhat similarly, an indicated raw or bulk mass flow rate MF ap parent (i.e., a mass flow 

30 rate of the entire two-phase flow) measured by a coriolis flowmeter may be different by a 

predictable or characterizable amount from an actual bulk mass flow rate MF^e. It should be 
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understood that correction techniques for corrected bulk mass flow rate MFtme may be different 
than the techniques for correcting for density. For example, various techniques for correcting a 
measured MF a p par ent to obtain an actual MFtme (or, at least, MF CO rrected) are discussed in U.S. Patent 
No. 6,505,519. 

5 Examples of detailed techniques for correcting p apP arent and MF ap p ar ent are discussed in 

more detail below. Generally speaking, though, with respect to FIG. 2, the digital transmitter is 
shown as including a density correction system 240, which has access to a density correction 
database 245, and a mass flow rate correction system 250, which has access to a mass flow 
correction database 255. As discussed in more detail below, the databases 245 and 255 may 

10 contain, for example, correction algorithms that have been derived theoretically or obtained 
empirically, and/or correction tables that provide corrected density or mass flow values for a 
given set of input parameters. The databases 245 and 255 also may store a variety of other types 
of information that may be useful in performing the density or mass flow corrections. For 
example, the density correction database may store a number of densities pii qu id corresponding to 

15 particular liquids (e.g., water or oil). 

Further in FIG. 2, a void fraction determination/correction system 260 is operable to 
determine a void fraction of a two-phase flow including a liquid and a gas. In one 
implementation, for example, the void fraction determination/correction system 260 may 
determine an actual void fraction a tru e from the corrected density p^. In another 

20 implementation, the void fraction determination/correction system 260 may input an apparent or 
indicated void fraction measurement obtained by the void fraction sensor 235, and may correct 
this measurement based on an error characterization similar to the density and mass flow 
techniques referred to above. In another implementation, the void fraction sensor 235 may be 
operable to directly measure an actual void fraction a^e, in which case the void fraction 

25 determination/correction system 260 simply inputs this measurement. 

Once the factors of p true , MFtme, and atme have been determined, and perhaps in 
conjunction with other known or discoverable quantities, a flow component mass flow rate 
determination system 265 operates to simultaneously determine a mass flow rate for the liquid 
phase component and a mass flow rate for the gas phase component. That is, the transmitter 104 

30 is operable to determine individual flowrates MFij qu jd and MF gas of the flow components, as 

opposed to merely determining the bulk flowrate of the combined or total two-phase flow MFtme- 
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Although, as just referred to, such measurements may be determined and/or output 
simultaneously, they also may be determined separately or independently of one another. 

Once the component flow rates MF|j qu jd and MF gas have been determined in the manner 
generally outlined above, these initial determinations may be improved upon by a process that 
5 relies on superficial velocities of the flow components, slip velocities between the components, 
and/or an identified flow regime of the flow. In this way, improved values for flow rates MFij qu i d 
and MFg as may be obtained, or may be obtained over time as those flow rates change. 

Superficial velocities are referred to herein as those velocities that would exist if the same 
mass flow rate of a given phase was traveling as a single phase through the flowtube 215. A 
10 superficial velocity determination/correction system 270 is included in the transmitter 104 for, 
for example, determining an apparent or corrected superficial velocity of a gas or liquid in the 
two-phase flow. 

Slip velocities refer to a condition in which gas and liquid phases in a two-phase flow 
have different average velocities. That is, an average velocity of a gas AV gas is different from an 
15 average velocity of a liquid AVij qu j d . As such, a phase slip S may be defined as S = 
AV gas /AVii quid . 

A flow regime is a term that refers to a characterization of the manner in which the two 
phases flow through the flowtube 215 with respect to one another and/or the flowtube 215, and 
may be expressed, at least partially, in terms of the superficial velocities just determined. For 

20 example, one flow regime is known as the "bubble regime," in which gas is entrained as bubbles 
within a liquid. As another example, the "slug regime" refers to a series of liquid "plugs" or 
"slugs" separated by relatively large gas pockets. For example, in vertical flow, the gas in a slug 
flow regime may occupy almost an entire cross-sectional area of the flowtube 215, so that the 
resulting flow alternates between high-liquid and high-gas composition. Other flow regimes are 

25 known to exist and to have certain defined characteristics, including, for example, the annular 
flow regime, the dispersed flow regime, and froth flow regime, and others. 

The existence of a particular flow regime is known to be influenced by a variety of 
factors, including, for example, a gas void fraction in the fluid flow, an orientation of the 
flowtube 215 (e.g., vertical or horizontal), a diameter of the flowtube 215, the materials included 

30 within the two-phase flow, and the velocities (and relative velocities) of the materials within the 
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two phase flow. Depending on these and other factors, a particular fluid flow may transition 

between several flow regimes over a given period of time. 

Information about phase slip may be determined at least in part from flow regime 

knowledge. For example, in the bubble flow regime, assuming the bubbles are uniformly 
5 distributed, there may be little relative motion between the phases. Where the bubbles 

congregate and combine to form a less uniform distribution of the gas phase, some slippage may 

occur between the phases, with the gas tending to cut through the liquid phase. 

In FIG. 2, a flow regime determination system 275 is included that has access to a 

database 280 of flow regime maps. In this way, information about an existing flow regime, 
10 including phase slip information, may be obtained, stored, and accessed for use in 

simultaneously determining liquid and gas mass flow rates within a two-phase flow. 

In FIG. 2, it should be understood that the various components of the digital transmitter 

104 are in communication with one another, although communication links are not explicitly 

illustrated, for the sake of clarity. Further, it should be understood that conventional components 
15 of the digital transmitter 104 are not illustrated in FIG. 2, but are assumed to exist within, or be 

accessible to, the digital transmitter 104. For example, the digital transmitter 104 will typically 

include (bulk) density and mass flow rate measurement systems, as well as drive circuitry for 

driving the driver 210. 

FIG. 3 is a flowchart 300 illustrating an operation of the coriolis flowmeter 200 of FIG. 2. 

20 Specifically, FIG. 3 illustrates techniques by which the flowmeter 200 of FIG. 2 is operable to 
simultaneously determine liquid and gas flow rates MFii quid and MF gas for a two-phase flow. 

In FIG. 3, it is determined that a gas/liquid two-phase flow exists in the flowtube 215 
(302). This can be done, for example, by an operator during configuration of the mass 
flowmeter/densitometer for gas/liquid flow. As another example, this determination may be 

25 made automatically by using a feature of the coriolis meter to detect that a condition of two- 
phase gas-liquid flow exists. In the latter case, such techniques are described in greater detail in, 
for example, U.S. Patent number 6,31 1,136 and U.S. Patent number 6,505,519, incorporated by 
reference above. 

Once the existence of two-phase flow is established, a corrected bulk density p^e is 
30 established (304) by the density correction system 240, using the density correction database 245 
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of the transmitter 104. That is, an indicated density p apP arent is corrected to obtain ptr Ue . 
Techniques for performing this correction are discussed in more detail below. 

Once ptrue is determined, a corrected gas void fraction atrue may be determined (306) by 
the void fraction determination/correction system 260. Also, a corrected bulk mass flow rate 
5 MFtrue is determined (308) by the mass flow rate correction system 250. As with density, 

techniques for obtaining the corrected void fraction a^e and mass flow rate MFtrue are discussed 
in more detail below. 

In FIG. 3, it should be understood from the flowchart 300 that the determinations of ptrue, 
atrue, and MFtrue may occur in a number of sequences. For example, in one implementation, the 

10 corrected void fraction a tru e is determined based on previously-calculated corrected density p^e, 
whereupon the corrected mass flow rate MF^e is determined based on a^ue- In another 
implementation, a^e and p trU e may be calculated independently of one another, and/or p^e and 
MFtme may be calculated independently of one another. 

Once corrected density p^e, corrected void fraction atrue, and corrected mass flow rate 

15 MRn-ue are known, then the mass flow rates of the gas and liquid components are determined 
(310) by the flow component mass flow rate determination system 265. Techniques for 
determining the liquid/gas component flow rates are discussed in more detail below with respect 
to FIG. 4. 

Once determined, the liquid/gas component flow rates may be output or displayed (312) 
20 for use by an operator of the flowmeter. In this way, the operator is provided, perhaps 

simultaneously, with information about both the liquid mass flow rate MFi iqui d and the gas mass 

flow rate MF gas of a two-phase flow. 

In some instances, this determination may be sufficient (314), in which case the 

outputting of the liquid/gas component flow rates completes the process flow. However, in other 
25 implementations, the determination of the individual component mass flow rates may be 

improved upon by factoring in information about, for example, the superficial velocities of the 

gas/liquid components, the flow regime(s) of the flow, and phase slip, if any, between the 

components. 

In particular, superficial velocities of the gas and liquid, SV gas and SVii quid are determined 
30 as follows. Gas superficial velocity SV gas is defined as: 
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SVgas = MF gas /(pgas * A T ) 



Eq. 1 



where the quantity At represents a cross-section area of the flowtube 215, which may be taken at 
a point where a void fraction of the flow is measured. Similarly, a liquid superficial velocity 
SVii qu id is defined as: 



As shown in Eqs. 1 and 2, determination of superficial velocities in this context relies on 
the earlier determination of MF gas and MFij qu jd. It should be understood from the above 
description and from FIG. 3 that MF gas and MFii qu id represent corrected or true mass flow rates, 
MF'™ and MF^ W since these factors are calculated based on ptrue, octme? and MF^e. As a result, 

the superficial velocities SV gas and SVi iquid represent corrected values SV™ and SVj™ rf . 

Further, the density values p gas and pij qu id refer, as above, to known densities of the liquid and gas 
in question, which may be stored in the density correction database 245. As discussed in more 
detail below with respect to techniques for calculating corrected density p trU e, the density values 
p gas and piiquid niay be known as a function of existing temperature or pressure, as detected by 
temperature sensor 220 and pressure sensor 225. 

Using the superficial velocities and other known or calculated factors, some of which 
may be stored in the flow regime maps database 280, a relevant flow regime and/or phase slip 
may be determined (318) by the flow regime determination/correction system 275. Once 
superficial velocities, flow regime, and phase slip are known, further corrections may be made to 
the corrected bulk density p tru e, corrected bulk mass flow rate MFtrue, and/or corrected void 
fraction ow In this way, as illustrated in FIG. 3, component flow rates MF gas and MFi iqu i d may 
be determined. 

Flow regime(s) in two phase liquid/ gas flow may be described by contours on a graph 
plotting the liquid superficial velocity versus the gas superficial velocity. As just described, an 
improvement to determinations of p tr ue> (Xtrue, and/or MF^e may be obtained by first establishing 
an approximate value of the liquid and gas flow rates, and then applying a more detailed model 
for the flow regime identified. For example, at relatively low GVF and relatively high flow there 



SVii qu id = MFii qui d/(pii q uid * A T ) 



Eq.2 
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exists a flow regime in which the aerated fluid behaves as a homogenous fluid with little or no 
errors in both density and mass flow. This can be detected as homogenous flow requiring no 
correction, simply using observation of the drive gain, which shows little or no increase in such a 
setting, despite a significant drop in observed density. 
5 FIG. 4 is a flowchart 400 illustrating techniques for determining liquid and gas flow rates 

MFii qu id and MF gas for a two-phase flow. That is, the flowchart 400 generally represents one 
example of techniques for determining liquid and gas flow rates (310), as described above with 
respect to FIG. 3. 

In FIG. 4, the determination of liquid and gas flow rates (310) begins with inputting the 
10 corrected density, void fraction, and mass flow rate factors pu-ue, ot tr ue> and MFtrue (402). In a first 
instance, (404), the liquid and gas flow rates are determined (406) using Eqs. 3 and 4: 

MFgas = Cttrue (pgas/ptrue) (MF true ) Eq. 3 

MFii quid = (I" a trU e)(pliquid/ptrue) (MF^e) Eq. 4 

15 

Eqs. 3 and 4 assume that there is no slip velocity (i.e., phase slip) between the liquid and 
gas phases (i.e., average velocity of the gas phase, AV gas , and average velocity of the liquid 
phase, AVii qu i d , are equal). This assumption is consistent with the fact that, in the first instance, 
superficial velocities and flow regimes (and therefore, phase slip) have not been determined. 

20 In the second instance and thereafter (404), a determination is made, perhaps by the flow 

regime determination/correction system 275, as to whether phase slip exists (408). If not, then 
Eqs. 3 and 4 are used again (406) or the process ends. 

If phase slip does exist (408), defined above as S = AV gas /AViiq U i d , the terms MF gas and 
MFii qu id are calculated using the cross-sectional area of the flowtube 215, A T , as also used in the 

25 calculation of superficial velocities in Eqs. 1 and 2 (410). Using the definition of slip S just 
given, 

MFgas = Pgas (Ottrue A T ) (AV gas) - p gas ((Xtrue AT)(S)(AV liqu i d ) Eq. 5 

MFiiq U i d = Pliquid ((1-CCtrue) A T ) (AVi iqu i d ) Eq. 6 

30 

Since MF^e = MF gas + MF]j qu id, Eqs. 5 and 6 may be solved for AVi iqu id to obtain Eq. 7: 
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AV liquid = MFtrue / (A T (pgasOttrueS + Piiquid(l-Ottrue))) 



Eq. 7 



As a result, the liquid and gas flow rates are determined (406) using Eqs. 8 and 9: 



MFii qu id = [piiquid (l-a tr ue)/( PgasOCtrueS + pliquid (1-OCtrue))] [MFtrue] 
MF gas = MFt ru e " MFi iquid 



Eq. 8 
Eq.9 



As described above, gas entrained in liquid forms a two-phase flow. Measurements of 
such a two-phase flow with a Coriolis flowmeter result in indicated parameters p apP arent> ot apP arent, 
and MFapparent for density, void fraction, and mass flow rate, respectively, of the two-phase flow. 
Due to the nature of the two-phase flow in relation to an operation of the Coriolis flowmeter, 
these indicated values are incorrect by a predictable factor. As a result, the indicated parameters 
may be corrected to obtain actual parameters pn-ue, ^true, and MF^e- In turn, the actual, corrected 
values may be used to simultaneously determine individual flow rates of the two (gas and liquid) 
components. 

FIGS. 5A and 5B are graphs illustrating a percent error in a measurement of void fraction 
and liquid fraction, respectively. In FIG. 5A, the percent error is a density percent error that is 
dependent on various design and operational parameters, and generally refers to the deviation of 
the apparent (indicated) density from the true combined density that would be expected given the 
percentage (%) of gas in liquid. 

In FIG. 5B, true liquid fraction versus indicated liquid fraction is illustrated. FIG. 5B 
shows the results, for the relevant flowmeter design, of several line sizes and flow rates. In 
more general terms, the functional relationship may be more complex and depend on both line 
size and flowrate. In FIG. 5B, a simple polynomial fit is shown that can be used to correct the 
apparent liquid fraction. 

Other graphing techniques may be used; for example, true void fraction may be plotted 
against indicated void fraction. For example, FIG. 6 is a graph illustrating a mass flow error as a 
function of a drop in density for a flowtube having a particular orientation and over a selected 
flow range. 
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FIG. 7 is a flowchart 700 illustrating techniques for correcting density measurements 
(304 in FIG. 3). In FIG. 7, the process begins with an inputting of the type of flowtube 215 
being used (702), which may include, for example, whether the flowtube 215 is bent or straight, 
as well as other relevant facts such as a size or orientation of the flowtube 215. 
5 Next, a gas-free density of the liquid, pii qu id is determined (704). This quantity may be 

useful in the following calculation(s), as well as in ensuring that that other factors that may 
influence the density measurement p apP arent, such as temperature, are not misinterpreted as void 
fraction effects. In one implementation, the user may enter the liquid density pii quid directly, 
along with a temperature dependence of the density. In another implementation, known fluids 
10 (and their temperature dependencies) may be stored in the density correction database 245, in 
which case the user may enter a fluid by name. In yet another implementation, the flowmeter 
200 may determine the liquid density during a time of single-phase, liquid flow, and store this 
value for future use. 

An indicated mass flow rate MF ap parent is read from the Coriolis meter (706), and then an 
15 indicated density p app arent is read from the Coriolis meter (708). Next, the density correction 

system 240 applies either a theoretical, algorithmic (710) or empirical, tabular correction (712) to 
determine the true density p trU e of the gas/liquid mixture. The quantity ptrue may then be output 
as the corrected density (714). 

An algorithmic density correction (710) may be determined based on the knowledge that, 
20 if there were no effect of the two-phase flow from the normal operation of a Coriolis meter when 
used to measure density, the indicated density would drop by an amount derived from the 
equation describing void fraction, which is set forth above in terms of volume flow and repeated 
here in terms of density as Eq. 10: 

25 Gt(%) = [(Pa PP arent - Piiquid) / ( Pgas ~ piiquid)] X 100 Eq. 10 

This can be used to define a quantity "density drop," or A p, as shown in Eq. 11: 

Ap = (papparent - Piiquid) = <*(%) X ( p gas - Piiquid) / 100 Eq. 1 1 

30 
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Note that Eq. 1 1 shows the quantity A p as being positive; however, this quantity could 
be shown as a negative drop simply by multiplying the right-hand side of the equation by -1, 
resulting in Eq. 12: 

5 Ap= ( pliquid - Papparent) = 0C(%) X ( Pliquid ~ Pgas ) / 100 Eq. 12 

The quantity p gas may be small compared to pii qu id> in which case Eq. 12 may be 
simplified to Eq. 13: 

10 A p= (Pliquid- Papparent ) = <*(%) * pi iquid / 100 Eq. 13 

As discussed extensively above, density measurements by a Coriolis meter, or any 
vibrating densitometer, generally are under-reported by the meter, and require correction. 
Accordingly, under two-phase flow Eqs. 12 or 13 may thus be used to define the following two 
15 quantities: a corrected or true density drop, A p^e, and an indicated or apparent density drop, 
A Papp. Using Eq. 13 as one example, this results in Eqs. 14 and 15: 

A ptrue = ( Pliquid^ Ptrue) = <*(%) X pii quid / 100 Eq. 14 

20 A pa P p = ( pliquid- Papparent ) = <*(%) X pi iquid / 100 Eq. 15 

There can be derived or empirically determined a relationship between A ptrue and A 
papparent and apparent mass flow rate, MF a pparent, as well as other parameters, such as, for example, 
drive gain, sensor balance, temperature, phase regime, etc.). This relationship can be expressed 
25 as shown as A ptrue = f (MF appa rent, drive gain, sensor balance, temperature, phase regime, and/or 
other factors). 

As a result, the relationship may generally be derived, or at least proven, for each 
flowtube in each setting. For one model flowtube, known and referred to herein as the 
Foxboro/Invensys CFS10 model flowtube, it has been empirically determined that for some 
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conditions the above functional relationship can be simplified to be only a function A p apP arent and 
of the form shown in Eq. 16: 



M 

^Ptrue = Z U i Apparent )' E q. 1 6 



i=0 



To force the condition for both sides of Eq. 16 to be zero when there is no apparent 
density drop relationship results in Eq. 17: 



M 

A Ptrue = Y, U i^P apparent Eq. 17 



M generally depends on the complexity of the empirical relationship, but in many cases can be as 
small as 2 (quadratic) or 3 (cubic). 

Once the true density drop is determined, then working back through the above equations 
it is straightforward to derive the true mixture density ptrue, as well as the true liquid and gas 
(void) fractions (the latter being discussed in more detail with respect to FIG. 9). 

A tabular correction for density (712) may be used when, for example, a functional 
relationship is too complex or inconvenient to implement. In such cases, knowledge of the 
quantities Apparent and AMF a pp arent may be used to determine Apt^e by employing a table having 
the form of a table 800 of FIG. 8. 

The table 800 may be, for example, a tabular look-up table that can be, for example, 
stored in the database 245, or in another memory, for use across multiple applications of the 
table. Additionally, the table may be populated during an initialization procedure, for storage in 
the database 245 for an individual application of the table. 

It should be understood that either or both of the algorithmic and tabular forms may be 
extended to include multiple dimensions, such as, for example, gain, temperature, balance, or 
flow regime. The algorithmic or tabular correction also may be extended to include other surface 
fitting techniques, such as, for example, neural net, radical basis functions, wavelet analyses, or 
principle component analysis. 
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As a result, it should be understood that such extensions may be implemented in the 
context of FIG. 3 during the approach described therein. For example, during a first instance, 
density may be determined as described above. Then, during a second instance, when a flow 
regime has been identified, the density may be further corrected using the flow regime 
information. 

FIG. 9 is a flowchart 900 illustrating techniques for determining void fraction 
measurements (306 in FIG. 3). In FIG. 9, the process begins with an inputting by the void 
fraction determination system 240 of the previously-determined liquid and bulk (corrected) 
densities, p| iquid and ptrue (902). 

A density of the gas, p gas is then determined (904). As with the liquid density pii qu id, there 
are several techniques for determining p gas . For example, p gas may simply be assumed to be a 
density of air, generally at a known pressure, or may be an actual known density of the particular 
gas in question. As another example, this known density p gas may be one of the above factors 
(i.e., known density of air or the specific gas) at an actual or calculated pressure, as detected by 
the pressure sensor 225, and/or at an actual or calculated temperature, as detected by the 
temperature sensor 220. The temperature and pressure may be monitored using external 
equipment, as shown in FIG. 2, including the temperature sensor 220 and/or the pressure sensor 
225. 

Further, the gas may be known to have specific characteristics with respect to factors 
including pressure, temperature, or compressibility. These characteristics may be entered along 
with an identification of the gas, and used in determining the current gas density p gas . As with 
the liquid(s), multiple gasses may be stored in memory, perhaps along with the characteristics 
just described, so that a user may access density characteristics of a particular gas simply by 
selecting the gas be name from a list. 

Once the factors pii qu id, Pgas, and p^e are known, then it should be clear from Eq. 10 that 
void fraction a true may be easily determined (906). Then, if needed, liquid fraction may be 
determined (908) simply by calculating l-atme- 

Although the above discussion presents techniques for determining void fraction atrue 

based on density, it should be understood that void fraction may be determined by other 

techniques. For example, an indicated void fraction a appar ent may be directly determined by the 

Coriolis flowmeter, perhaps in conjunction with other void fraction determination systems 
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(represented by the void fraction sensor 235 of FIG. 2), and then corrected based on empirical or 
derived equations to obtain atme- In other implementations, such external void fraction 
determining systems may be used to provide a direct measurement of a tr ue. 

FIG. 10 is a flowchart 1000 illustrating techniques for determining corrected mass flow 
5 rate measurements (308 in FIG. 3). In FIG. 10, the mass flow rate correction system 250 first 

inputs the previously-calculated corrected density drop Apt™ (1002), and then inputs a measured, 
apparent mass flow rate MF a p par ent (1004). 

The mass flow rate correction system 250 applies either a tabular (1006) or algorithmic 
correction (1008) to determine the true mass flow rate MFtrue of the gas/liquid mixture. The 
10 quantity MFtrue may then be output as the corrected mass flow rate (1010). 

In applying the tabular correction for mass flow rate (1006), knowledge of the quantities 
Aptrue and AMF a pp a rent may be used to determine MFtrue by employing a table having the form of a 
table 11 00 of FIG. 11. 

The table 1 100, as with the table 800 may be, for example, a tabular look-up table that 
15 can be, for example, stored in the database 245, or in another memory, for use across multiple 
applications of the table. Additionally, the table may be populated during an initialization 
procedure, for storage in the database 255 for an individual application of the table. 

Normalized values MF nor m_app and MF norm _true niay be used in place of the actual ones 
shown above, in order to cover more than one size coriolis flowtube. Also, the entries can be in 
20 terms of the correction, where the correction is defined by Eq. 18: 

AMF = MF^e - MFapparent Eq. 1 8 

The values in Eq. 18 should be understood to represent either actual or normalized values. 
25 In an algorithmic approach, as with density, the correction for mass flow may be 

implemented by way of a theoretical or an empirical functional relationship that is generally 

understood to be of the form AMF = f (MFapparent, void fraction, drive gain, sensor balance, 

temperature, phase regime, and/or other factors). 

For some cases the function can simplify to a polynomial, such as, for example, the 
30 polynomial shown in Eq. 19: 
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AMF = J ia ( .^(A P L)(MF^_ w ) Eq. 19 

1-0 7=0 

For some set of conditions, the functional relationship can be a combination of a 
polynomial and exponential, as shown in Eq. 20: 

AMF = aide (^ 2 ^a 4m ^a 5 rn) + ^2 + ^ + ^2 + ^ 

Eq. 20 

In Eq. 20, d = A p tru e, and m = f (MF ap parent). 

In one implementation, m in Eq. 20 may be replaced by apparent superficial liquid 
velocity SVi iquid which is given as described above by Eq. 2 as SVn qu i d = MFi iqu i d /(piiquid * At). In 
this case, pii qu id and flowtube cross-section A T are known or entered parameters, and may be real- 
time corrected for temperature using, for example, the on-board temperature measurement device 
220 of the digital controller/transmitter 104. 

It should be understood that, as with the density corrections discussed above, either or 
both of the algorithmic and tabular forms may be extended to include multiple dimensions, such 
as, for example, gain, temperature, balance, or flow regime. The algorithmic or tabular 
correction also may be extended to include other surface fitting techniques, such as, for example, 
neural net, radical basis functions, wavelet analyses, or principle component analysis. 

As a result, it should be understood that such extensions may be implemented in the 
context of FIG. 3 during the approach described therein. For example, during a first instance, 
mass flow rate may be determined as described above. Then, during a second instance, when a 
flow regime has been identified, the mass flow rate may be further corrected using the flow 
regime information. 

All of the above functional relationships for mass flow rate may be restated using gas 
fraction (a) or liquid fraction (100- a) instead of density drop, as reflected in the table 1 100 of 
FIG. 11. Also, although the above described methods are dependent on knowledge of the 
corrected density drop A ptrue> it should be understood that other techniques may be used to 
correct an indicated mass flow rate. For example, various techniques for correcting mass flow 
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rate measurements of a two-phase flow are discussed in U.S. Patent 6,505,519, incorporated by 
reference above. 

Having described density, void fraction, and mass flow rate corrections above in general 
terms, for the purpose of, for example, simultaneously calculating individual flow component 
(phases) flow rates in a two-phase flow, the below discussion and corresponding figures provide 
specific examples of implementations of these techniques. 

FIGS. 12-14 are graphs illustrating examples of density corrections for a number of 
flowtubes. In particular, the examples are based on data obtained from three vertical water 
flowtubes, the flowtubes being: Vi \ 3 / 4 ", and 1" in diameter. 

More specifically, the 14" data was taken with a 0.15kg/s flow rate and a 0.30kg/s flow 
rate; the 3 / 4 " data was taken with a 0.50kg/s flow rate and a 1 .00kg/s flow rate; and the 1" data 
was taken with a 0.50kg/s flow rate, a 0.90kg/s flow rate, and a 1.20kg/s flow rate. FIG. 12 
illustrates an error, e d , of the apparent density of the fluid-gas mixture (two-phase flow) versus 
the true drop in density of the fluid-gas mixture, A p^: 



A Anje =100- A '^ P,rue Eq.21 

P liquid 

ed=m . p w°™<- p <™ Eq 22 

Ptrue 



where, as above, p i iquid is the density of the gas-free liquid, p^e is the true density of the liquid- 
gas mixture, and p ap parent is the apparent or indicated density of the liquid-gas mixture. 

In FIG. 12, the correction is performed in terms of the apparent drop in mixture density, 
A Papparent, as shown in Eq. 23: 



^Papparent = H>0 ■ ^ Eq. 23 



apparent 



In FIG. 12, when fitting the data, both the apparent and true drop in density of the 
mixture were normalized to values between 0 and 1 by dividing them through by 100, where this 
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normalization is designed to ensure numerical stability of the optimization algorithm. In other 
words, the normalized apparent and true drop in mixture density are the apparent and true drop in 
mixture density defined as a ratio, rather than as a percentage, of the liquid density pii qu id> as 
shown in Eq. 24: 



k normalized ^P apparent p 
apparent ~~ \QQ 



The model formula, based on Eq. 17, provides Eq. 25: 



A ^normalized _ / a ^normalized \ 3 , /a ^normalized \ 2 . „ / a ^normalized \ T7- 
A A™ = a A^P apparent ) + <*2 ( A /W„, ) + ^3 i^P apparent ) Eq« 25 



In this case, the coefficients are a;= -0.51097664273685, 1.26939674868129, and 03= 

0.24072693 1 19420. FIGS. 13 A and 13B illustrate the model with the experimental data and 
the residual errors, as shown. FIGS. 14A and 14B give the same information, but with each flow 
rate plotted separately. 

To summarize, the drop in density correction is performed in the transmitter 104 by 
calculating the apparent density drop A papparem, using the apparent density value p apP arent and the 
liquid density pijquid- The value of the apparent drop in density is normalized to obtain 

^PaZarlnt^ = ^° p ^ ent ' so as ex pl a i ne d above, the drop in density is calculated as a ratio 

rather than a percentage. The density correction model(s) may then be applied to obtain the 
normalized corrected drop in mixture density Ap"^ r e mahzed . Finally, this value is un-normalized to 

obtain the corrected drop in density Ap^ = 100 • Ap™ r e mall2ed . Of course, the final calculation is 

not necessary if the corrected drop in mixture density Ap^ is defined as a ratio rather than 

percentage of the true value. 

FIGS. 15-20 are graphs illustrating examples of mass flow rate corrections for a number 
of flowtubes. In particular, the examples are based on data obtained from three vertical water 
flowtubes, the flowtubes being: J4", 3 / 4 ", and 1" in diameter. More specifically, the 54" data was 
taken with a 0.15kg/s flow rate and a 0.30kg/s flow rate; the 3 / 4 " data was taken with a 0.50kg/s 
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flow rate and a l.OOkg/s flow rate; and the 1" data was taken with 18 flow rates between 0.30kg/s 
and 3.0kg/s flow rate, with a maximum drop in density of approximately 30%. 

FIGS. 15A and 15B illustrate apparent mass flow errors for the data used to fit the model 
versus corrected drop in mixture density Ap^ and normalized true superficial fluid velocity; 
i.e., the apparent mass flow error curves per flowline, together with a scatter plot of the apparent 
mass flow error versus corrected drop in density Ap true and normalized true superficial fluid 
velocity v m , as shown in Eq. 26: 

v,„=^-,v,= M ' Eq.26 

V max P liquid * A T 

where m t is the true fluid mass flow, i.e. the value of the mass flow independently measured, 

p ii qu id is the liquid density, A T is the flowtube cross-section area, and v max is the maximum 

value for the superficial fluid velocity (here considered 12), so that v tn gives the ratio of the true 

superficial fluid velocity from the whole range of the flowtube 215. In these examples, both 
drop in mixture density and superficial fluid velocity are normalized between 0 and 1 prior to 
fitting the model, for the purpose of ensuring numerical stability for the model optimization 
algorithm. 

FIG. 16 illustrates apparent mass flow errors versus corrected drop in mixture density and 
normalized apparent superficial fluid velocity, with safety bounds for the correction mode. That 
is, FIG. 16 gives the scatter plot of the apparent mass flow errors versus corrected drop in density 

V 171 

and, this time, normalized apparent superficial fluid velocity v n = = , where m is 

v v • o - A 

max max H yA 

the apparent fluid mass flow (i.e. as measured by the transmitter 104). Superimposed on the plot 
are the boundaries defining the safe region for the model, i.e., the region for which the model is 
expected to give an accuracy similar with the one for the fit data. Using this nomenclature, the 

apparent mass flow error e is given by e = 100 • m 



The model formula for this situation is shown as Eq. 27: 

e n = a x dd cn • + a ^ + + fl|V a + ^ Eq . 27 
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where 

e m-m. ^ ^ 

100 m t 

where, in Eqs. 27 and 28, dd cn is the normalized corrected drop in mixture density, and 
v n is the normalized apparent superficial velocity of the liquid. 
5 In this case, the coefficient are: aj = -4.78998578570465, a2= 4.20395000016874, a 3 

= -5.93683498873342, a 4 = 12.03484566235777, = -7.70049487145105,^ = 
0.69537907794202, 07 = -0.52153213037389, ag= 0.36423791515369, mda 9 = - 
0.16674339233364 

FIG. 17 illustrates a scatter plot for the model residuals, together with the model formula 
10 and coefficients; i.e., shows model residuals versus the corrected drop in mixture density and 
normalized true fluid velocity. FIGS. 18A-18D and Figures 19A-19D give the model residual 
errors for the whole data set used to fit the model and the actual data alone, respectively. Finally, 
FIGS. 20A and 20B illustrate the model surface both interpolating and extrapolating outside the 
safe fit area. From FIGS. 16, 20A, and 20B, the apparent mass flow (superficial liquid velocity) 
15 and drop in density bounds for the model should be understood. 

To summarize, mass flow correction in the transmitter 104 is undertaken in this example 
by calculating an apparent drop in density, correcting it using the method(s) described above, and 
normalizing the resulting value by dividing it by 100 (or use the obtained normalized corrected 
drop in density from the density model). Then, a normalized superficial fluid velocity v n is 
20 calculated, and the model is applied to obtain an estimation of the normalized mass flow error 
e n , where this value gives the error of the apparent mass flow as a ratio of the true mass flow. 
The obtained value may be un-normalized by multiplying it by 100, to thereby obtain the mass 
flow error as a percentage of the true mass flow. Finally, the apparent mass flow may be 

corrected with the un-normalized mass flow error /w. = . 

25 As will be appreciated, the above description has a wide range of applications to improve 

the measurement and correction accuracy of a coriolis meter during two phase flow conditions. 
In particular, the techniques described above are particularly useful in measurement applications 
where the mass flow of the liquid phase and the mass flow of the gas phase must be measured 
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and/or corrected to a high level of accuracy. One exemplary application is the measurement of 
the mass flow of the liquid phase and the measurement of the gas phase in oil and gas production 
environments. 

The above discussion is provided in the context of the digital flowmeter of FIG. 2. 
However, it should be understood that any vibrating or oscillating densitometer or flowmeter, 
analog or digital, that is capable of measuring multi-phase flow that includes a gas phase of a 
certain percentage may be used. That is, some flowmeters are only capable of measuring process 
fluids that include a gas phase when that gas phase is limited to a small percentage of the overall 
process fluid, such as, for example, less than 5%. Other flowmeters, such as the digital 
flowmeter(s) referenced above, are capable of operation even when the gas void fraction reaches 
40% or more. 

Many of the above-given equations and calculations are described in terms of density, 
mass flow rate, and/or void fraction. However, it should be understood that the same or similar 
results may be reached using variations of these parameters. For example, instead of mass flow, 
a volumetric flow may be used. Additionally, instead of void fraction, liquid fraction may be 
used. 

A number of implementations have been described. Nevertheless, it will be understood 
that various modifications may be made. Accordingly, other implementations are within the 
scope of the following claims. 
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